%1 数字基带信号
f=10;
x=sin(2*pi*f*t);

figure;
subplot(921)
plot(t,x);
xlabel('时间 (s)');
ylabel('幅度');
title('数字基带信号的时域谱');

X = fft(x); % 对时域信号进行傅立叶变换来获得频域表示。

% 计算对应的频率轴
fs = 1/(t(2)-t(1)); % 采样频率
N = length(x); % 信号长度
frequencies = (0:N-1) * fs / N; % 计算对应的频率

subplot(922)
plot(frequencies, abs(X)/N); % 使用stem函数绘制幅度谱
xlabel('频率 (Hz)');
ylabel('幅度');
title('数字基带信号的频域谱');

%2 载波
fc=100;
carrier=cos(2*pi*fc*t);
subplot(923);
plot(t,carrier);
xlabel('时间 (s)');
ylabel('幅度');
title('载波信号的时域谱');
Carrier = fft(carrier); % 对时域信号进行傅立叶变换来获得频域表示。

% 计算对应的频率轴
fs = 1/(t(2)-t(1)); % 采样频率
N = length(carrier); % 信号长度
frequencies = (0:N-1) * fs / N; % 计算对应的频率

subplot(924)
plot(frequencies, abs(Carrier)/N); % 使用stem函数绘制幅度谱
xlabel('频率 (Hz)');
ylabel('幅度');
title('载波信号的频域谱');

%3 数字调制
modulated_signal=x.*carrier;
subplot(925)
plot(t,modulated_signal);
xlabel('时间 (s)');
ylabel('幅度');
title('数字调制后的信号的时域谱');

Modulated_signal=fft(modulated_signal);

% 计算对应的频率轴
fs = 1/(t(2)-t(1)); % 采样频率
N = length(modulated_signal); % 信号长度
frequencies = (0:N-1) * fs / N; % 计算对应的频率

subplot(926)
plot(frequencies, abs(Modulated_signal)/N); % 使用stem函数绘制幅度谱
xlabel('频率 (Hz)');
ylabel('幅度');
title('数字调制后的信号的频域谱');

% 对已调信号进行解调
demodulated_signal = modulated_signal .* carrier;
subplot(927);
plot(t, demodulated_signal);
xlabel('时间 (s)');
ylabel('幅度');
title('解调后的信号的时域谱');
Demodulated_signal = fftshift(fft(demodulated_signal)/N); % 频域谱
subplot(9,2,8);
plot(frequencies, abs(Demodulated_signal));
ADemodulated_signal = (fft(demodulated_signal)/N); % 频域谱
subplot(9,2,9);
plot(frequencies, abs(ADemodulated_signal));




